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ABSTRACT  Time-resolved  magneto-optical  imaging  measure¬ 
ments  show  that  an  ac  current  enables  the  vortex  matter  in 
YBa2Cu307_x  thin  films  to  reorganize  into  two  coexisting 
steady  states  of  driven  vortex  motion  with  different  characteris¬ 
tics:  a  quasi- static  disordered  glassy  state  in  the  sample  interior 
and  a  dynamic  state  of  plastic  motion  near  the  edges.  Finite- 
element  calculations  consistent  with  the  critical  state  model 
show  good  agreement  with  the  measured  field  profiles  in  the 
quasi- static  state  but  predict  a  larger  hysteretic  behavior  in  the 
dynamic  state. 

PACS  74.25.Qt;  74.25.Ha;  74.25.Sv;  74.78.Bz 


1  Introduction 

Type-II  and  high-temperature  superconductors 
(HTS)  subjected  to  quasi-static  changes  of  both  an  applied 
magnetic  field  and  transport  current  have  long  been  studied 
using  different  experimental  methods.  Hall  sensor  array  [1], 
scanning  Hall  probe  [2]  and  magneto-optical  imaging  [3] 
measurements  have  revealed  interesting  features  of  the  mag¬ 
netic  vortex  penetration  in  superconducting  thin  films,  such 
as  hysteretic,  inhomogeneous  and  disordered  states  induced 
by  the  applied  field  and/or  the  transport  current.  These  phe¬ 
nomena  are  generally  understood  in  terms  of  the  critical  state 
theory  [4,5],  which  describes  the  behavior  of  the  magnetic 
vortices  in  quasi-static  changes  of  an  applied  magnetic  field  or 
a  transport  current.  However,  recent  experiments  have  shown 
that  the  vortex  mobility  can  be  strongly  affected  by  a  dynamic 
magnetic  field  or  ac  current  producing  a  number  of  novel 
effects,  including:  low  frequency  noise  [6, 7],  slow  voltage  os¬ 
cillations  [8],  a  history  dependent  dynamic  response  [9, 10] 
and  a  memory  effect  related  to  the  direction,  amplitude,  dura¬ 
tion,  and  frequency  of  the  applied  current  [11].  Some  of  these 
effects  in  vortex  physics  have  been  ascribed  to  dynamic  in¬ 
stabilities  [1, 10-12]  and  transient  [13, 14]  or  steady  states  of 
driven  vortex  motion  created  during  the  reordering  induced  by 
the  applied  current.  Even  more  interesting  is  the  dynamic  in- 
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teraction  of  the  vortex  matter  (VM)  with  a  magnetic  field  and 
a  current  both  applied  at  the  same  time,  which  is  relevant  to 
most  HTS  applications  [15]. 

Here  we  present  a  detailed  study  of  the  effects  of  an  ac 
current  and  a  magnetic  field  on  the  dynamical  vortex  states 
in  a  YBa2Cu307_x  (YBCO)  thin  film  by  using  time-resolved 
magneto-optical  (TRMO)  imaging.  Since  TRMO  imaging  of¬ 
fers  simultaneously  high  spatial  and  temporal  resolution,  it 
is  suitable  for  accurate  dynamic  studies  of  magnetic  flux  dis¬ 
tributions.  We  compare  cross-sectional  field  and  current  pro¬ 
files  from  our  measurements  with  simulations  using  a  recently 
developed  method  based  on  finite-element  (FEM)  calcula¬ 
tions  [16],  which  accounts  for  the  thermally  activated  flux 
creep.  Our  results  show  that  a  quasi-static  disordered  glassy 
state  develops  in  the  sample  interior  in  agreement  with  FEM 
simulations.  A  dynamic  state  of  plastic  motion  forms  near  the 
edges,  which  exhibits  negligible  hysteretic  behavior  and  re¬ 
duced  pinning  strength.  This  effect  is  not  revealed  in  the  FEM 
calculations,  which  rely  on  the  critical  state  model. 

2  Experimental 

Two-dimensional  maps  of  the  magnetic  field  per¬ 
pendicular  to  the  YBCO  thin  films  are  recorded  with  TRMO 
imaging.  A  detailed  description  of  the  setup  can  be  found 
in  [17].  Time  resolution  is  achieved  by  stroboscopic  imaging 
using  a  Q-switched  frequency-doubled  Nd:YLF  laser  (100-ns 
pulse  length,  527-nm  wavelength)  with  a  pulse  repetition  fre¬ 
quency  (PRF)  tunable  from  single  shot  to  1200  Hz.  The  PRF 
of  the  laser  is  synchronized  to  the  frequency  of  an  ac  current 
source.  By  shifting  the  phase,  TRMO  images  can  be  taken  at 
specific  phase  points  over  the  whole  period  of  the  ac  current 
with  an  accuracy  <  1  |xs.  We  use  a  polarizing  microscope  in 
a  cross-polarized  configuration,  a  digital  camera  with  a  reso¬ 
lution  of  1344  x  1024  pixels,  and  a  Bi-substituted  ferrite  gar¬ 
net  crystal  as  a  magneto-optical  (MO)  indicator.  In  the  TRMO 
measurements,  a  4x-objective  is  selected  to  obtain  a  wide  field 
of  view  of  3 .4  x  2.6  mm2,  which  restricts  the  spatial  resolution 
to  a  few  microns. 

The  YBCO  samples  were  grown  by  pulsed  laser  depo¬ 
sition  on  a  LaA103  or  SrTi03  substrate  with  typical  dimen¬ 
sions  of  10  mm  x  5  mm  x  250  nm  [18].  The  critical  transi¬ 
tion  temperature  ( Tc )  of  the  samples  is  91  K,  as  determined 
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FIGURE  1  Time  resolved  magneto-optical  im¬ 
ages  for  different  values  of  phase  (p  of  the  applied 
ac  current.  The  white  arrow  in  the  center  of  the 
images  indicates  the  current  direction.  The  ap¬ 
plied  magnetic  field  is  normal  to  the  surface.  The 
red  arrows  point  towards  the  flux/anti-flux  anni¬ 
hilation  front 


by  ac  susceptibility  measurements.  The  samples  are  bridged 
using  a  photolithographic  technique  with  a  length  of  6  mm 
and  a  cross-section  2 w  =  460  pun  to  reduce  their  critical  cur¬ 
rent.  A  nanoparticulate  dispersion  provides  strong  pinning 
centers  in  the  sample  yielding  a  critical  current  density  Jc  = 
3.4  x  106  A/cm2  at  77  K. 

In  the  TRMO  experiment,  the  sample  is  cooled  down  to 
T  =  24  K  in  a  zero  magnetic  field.  Upon  reaching  the  tem¬ 
perature,  a  magnetic  field  Fa  =  10  mT  is  applied  perpendicu¬ 
lar  to  the  sample  surface.  An  ac  current  I(t)  =  /o  s\n(2n ft) 
with  /o  =  8.54  A  and  /  =  1  kHz  is  sent  through  the  sample.  At 
this  temperature  the  applied  current  corresponds  to  about  32% 
of  the  critical  current.  We  collect  25  TRMO  images  over  the 
whole  period  of  the  ac  current  with  a  time  step  of  40  |xs. 

3  Results  and  discussion 

Figure  1  shows  six  images  with  focus  on  a  small  re¬ 
gion  of  the  bridge  of  about  10%  of  its  length,  which  are  taken 
at  different  phase  points  <p  =  2nft  during  the  whole  cycle. 
Bright  areas  in  the  images  correspond  to  regions  of  high  mag¬ 
netic  field  near  the  edges  of  the  sample.  The  dark  area  in  the 
center  of  the  images  indicates  the  flux-free  region.  The  uni¬ 
formly  gray  region  represents  the  intermediate  magnetic  field 
outside  the  sample.  The  penetration  of  the  magnetic  flux  lines 
in  the  sample  is  cyclic,  asymmetric  and  inhomogeneous,  es¬ 
pecially  near  the  edges  of  the  sample  where  the  vortices  form 
a  pattern  of  narrow  parallel  bands,  indicated  by  the  red  arrows 
in  Fig.  1 .  The  complex  distribution  of  vortices  specifically  vis¬ 
ible  in  the  images  at  (p  =  |  and  <p  =  ^  is  a  clear  indication  that 
the  ac  current  rearranges  the  VM  during  the  cycle. 

Figure  2  shows  cross-sectional  magnetic  field  and  current 
density  profiles  averaged  over  a  sample  length  of  about  1  mm 
for  five  selected  phase  points  between  <p  =  |  and  (p  =  ^  . 
Only  the  central  part  of  the  profiles  ( |  ^  |  =  \x\  <  2)  are  shown 
in  Fig.  2,  but  the  measured  profiles  extend  up  to  \x  \  = 5  out¬ 
side  the  sample.  Current  density  profiles  (Fig.  2b)  are  calcu¬ 
lated  from  the  magnetic  field  profiles  by  numerical  inversion 
of  the  Biot-Savart  law  [19].  The  gray  shaded  regions  in  Fig.  2 
indicate  the  dynamic  vortex  state. 

Figures  1  and  2  indicate  clearly  that  the  region  of  the  high¬ 
est  magnetic  field  changes  from  one  edge  (x  =  1)  to  the  other 
(x  =  —  1)  according  to  the  direction  of  the  applied  ac  current 
and  the  induced  self-field,  since  the  external  magnetic  field 
is  constant  during  the  cycle.  The  regions  where  vortices  and 
anti- vortices  annihilate  appear  dark  in  the  TRMO  images.  For 


a  x 


FIGURE  2  Field  (a)  and  current  density  (b)  profiles  for  different  values  of 
the  phase.  The  phase  points  are  indicated  by  the  inset.  The  dotted  lines  show 
the  sample  edges 


example  in  Fig.  1  at  <p  =  |  and  cp=^-,  two  such  regions  are 
visible  as  narrow  dark  lines  (marked  by  red  arrows).  These 
two  lines  constitute  the  boundary  of  the  negative  flux  region, 
and  we  accounted  for  them  during  the  calibration  of  the  im¬ 
ages  [19]. 

When  the  magnetic  field  Ba(z )  >  Hc\  is  applied  to  the 
superconducting  film,  vortices  penetrate  the  region  (\x\  ^  xp), 
where  xp  =  0.65  (vp  <  w)  is  the  field  dependent  penetration 
depth  as  indicated  in  Fig.  2.  This  induces  a  shielding  current 
density  Jsc  (x,  y )  circulating  in  the  sample  that  according  to  the 
critical  state  model  reaches  its  maximum  value  Jc  in  the  re¬ 
gion  (\x\  ^  vp)  where  the  flux  penetrates.  The  Lorentz  force, 
Fl  =  (JxB),  created  by  the  shielding  current,  is  responsible 
for  moving  the  vortices.  The  final  equilibrium  of  the  vortex 
phase  in  the  HTS  thin  film  is  determined  by  the  competition  of 
Fl  with  four  major  interactions  [20]: 
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-  The  thermal  fluctuations  (Ft)  that  favor  a  vortex  liquid  of 
lines  or  pancakes 

-  The  vortex-vortex  interaction  (Fvv)  that  favors  a  perfect 
ordered  vortex  lattice 

-  The  pinning  (Fp)  that  instead  favors  an  amorphous  or 
glassy  solid 

-  The  coupling  force  (Fc)  between  layers  that  forms  vortex 
lines  from  weakly  interacting  pancake  vortices  in  adjacent 
layers 

In  HTS  all  these  interactions  have  comparable  energy  scales 
so  that  the  VM  at  equilibrium  is  rich  in  different  phases.  In 
YBCO  thin  films  at  low  temperature  and  low  field  the  pin¬ 
ning  prevails  over  the  other  interactions,  due  to  the  strong 
anisotropy,  the  random  distribution  of  surface  barriers  and  de¬ 
fects.  The  VM  at  equilibrium  is  expected  to  be  disordered  and 
glassy;  however,  the  application  of  a  transport  current  density 
y)  at  a  fixed  temperature  can  rearrange  the  distribution 
of  vortices. 

/tr  increases  the  total  current  density  J  and  hence  FL  on 
one  side  of  the  sample  while  it  decreases  both  on  the  other 
side.  Additional  flux  penetrating  from  the  edge  shifts  the  flux 
penetration  front  further  into  the  sample  towards  the  center 
until  it  reaches  a  new  equilibrium  point  where  JA  Ft  =  0.  The 
critical  state  model  predicts  that  the  maximum  current  dens¬ 
ity  region  is  extended  uniformly  to  the  new  position  of  the  flux 
penetration  front  and  allows  the  calculation  of  the  new  cur¬ 
rent  and  field  distributions  [19].  However  rapid  variations  of 
the  field,  as  in  the  case  of  high-frequency  alternating  trans¬ 
port  currents,  can  lead  to  a  complete  rearrangement  of  the  VM 
and  a  more  complicated  behavior  than  predicted  by  the  critical 
state  model.  Metastable,  disordered  vortex  phases  created  by 
the  transport  current  or  a  pulsed  magnetic  field  have  been  pre¬ 
viously  observed  in  both  type-II  superconductors  [1,21]  and 
HTS  [11, 22].  In  this  case  analytical  solutions  describing  the 
time-evolution  of  the  system  become  challenging  to  compute 
or  even  impossible  to  obtain  especially  if  other  effects,  such 
as  vortex  phase  transitions  or  flux  instabilities,  are  included  in 
the  theoretical  description. 

In  order  to  study  the  effects  of  the  ac  current  on  the  VM 
we  simulated  the  phase  evolution  of  the  field  profiles  by  adopt¬ 
ing  a  recently  developed  method  based  on  FEM  calculations. 
The  model  uses  the  direct  magnetic  field  formulation  with¬ 
out  the  use  of  vector  or  scalar  potentials  (common  in  con¬ 
ventional  formulations)  and  relies  on  first-order  edge  finite 
elements  [16].  The  governing  equations  are  Maxwell  equa¬ 
tions  and  the  superconducting  material  is  described  by  means 
of  a  non-linear  E-J  characteristic  which  takes  into  account 
the  transition  from  the  superconducting  to  normal  state.  In 
particular,  the  resistivity  of  the  superconductor  is  expressed 
by 


(1) 


where  Ec  is  the  critical  electrical  field,  Jc  is  the  critical  cur¬ 
rent  density  and  n  is  the  power  index  defining  the  steepness  of 
the  transition.  In  our  simulations,  we  used  the  values  for  the 
YBCO  thin  films  used  in  the  experiment  at  T  =  24  K:  Jc  = 
2.2  x  1011  Am-2,  Ec  =  10-4  Vm-1  and  n  =  19.  The  non¬ 
linear  dependence  of  the  current  and  electric  field  expressed 


in  (1)  is  determined  by  the  thermally  activated  flux  creep  [20]. 
The  value  of  the  exponent  n  is  related  to  the  ratio  between 
the  thermal  energy  and  the  flux  pinning  potential  Uo,  which 
at  the  measured  temperature  corresponds  to  40  meV.  The  field 
and  the  ac  current  simulation  parameters  are  fixed  during  the 
simulation  to  experimental  values:  Ba  =  lOmT  for  the  mag¬ 
netic  field  and  /o  =  8.54  A  and  /  =  1  kHz  for  the  current.  In 
order  to  avoid  transient  effects,  the  tenth  cycle  of  simulation 
has  been  considered  for  analysis.  In  Fig.  2  we  compare  the 
calculated  profiles  for  two  phase  points  cp  =  |  and  cp  =  ^ 
to  the  measured  data.  The  FEM  simulations  reproduce  the 
flux  penetration  front,  the  shifting  of  the  maxima  with  the 
phase  and  the  asymptotic  behavior  outside  the  sample  for 
all  the  measured  phase  points.  Differences  can  be  observed 
near  the  edges  of  the  sample  where  the  maximum  peak  of 
the  measured  profiles  is  reduced  in  intensity  and  the  peak 
at  the  opposite  edge  is  shifted  further  into  the  center  of  the 
sample. 

Figure  3  shows  the  field  (Fig.  3a)  and  current  (Fig.  3b)  dis¬ 
tribution  profiles  measured  and  simulated  with  no  net  trans¬ 
port  current  flowing  in  the  sample.  The  curves  indicated  by 
blue  squares  and  red  triangles  are  the  profiles  measured  at 
(p  =  0  and  cp  =  7i  when  the  voltage  during  the  ac  cycle  is 
zero.  At  these  phase  points  the  electrons  flow  in  response  to 
the  trapped  magnetic  flux  lines  and  the  measured  current  ef¬ 
fectively  represents  a  shielding  current.  These  two  current 


FIGURE  3  Current  density  profiles  corresponding  to  the  shielding  current 
at  <p  =  0  and  <p  =  n,  compared  to  the  shielding  current  density  for  the  case  of 
only  an  applied  magnetic  held  of  10  mT  without  transport  current.  The  dotted 
lines  indicate  the  sample  edges 
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profiles  overlap  fairly  well  over  the  whole  sample  width, 
with  the  exception  of  a  small  deviation  in  a  region  close 
to  the  edges.  For  comparison,  we  also  show  the  shielding 
current  distribution  (black  circle)  induced  by  only  the  ap¬ 
plied  field  #a  =  10  mT.  The  FEM  simulations  are  also  plot¬ 
ted  in  Fig.  3  as  blue  and  red  solid  lines  corresponding  to 
cp  =  0  and  <p  =  7T,  respectively.  The  simulations  at  these  phase 
points  show  a  stronger  hysteretic  behavior  than  the  meas¬ 
ured  profiles  near  the  sample  edge.  The  effect  of  the  ac 
transport  current  is:  (1)  to  reduce  the  peak  shielding  current, 
(2)  to  move  the  shielding  current  further  inside  the  sample, 
and  (3)  to  create  a  narrow  gap  at  the  edges  where  /sc(^  = 
0, 7r)  is  strongly  reduced  and  negligible  hysteretic  effects  are 
observed. 

The  applied  current  softens  the  pinning,  allowing  the  vor¬ 
tices  to  find  more  strongly  pinned  states  inside  the  sample. 
The  current  that  flows  at  the  edges  injects  flux  creating  a  tran¬ 
sient  disordered  state.  After  a  few  cycles,  subsequent  to  the 
flux  injection  the  VM  reorganizes  in  a  steady  state  of  driven 
vortex  motion.  The  motion  occurs  in  the  presence  of  strong 
pinning  centers  that  impede  “smooth”  entrance  of  the  injected 
fluxions.  At  low  temperatures  there  is  a  stationary  phase  at 
a  low  driving  force  and  a  dynamic  phase  above  a  critical 
de-pinning  force.  Since  the  penetration  and  exit  of  vortices 
requires  a  large  FL,  due  to  the  surface  barriers  much  of  the 
applied  transport  current  flows  at  the  edges  in  order  to  pro¬ 
vide  the  necessary  driving  force  [20, 23]  and  the  rest  phase  is 
established  in  the  center  of  the  sample.  The  nature  of  the  mo¬ 
tion  above  the  critical  force  can  exhibit  a  wide  variation  [24]. 
In  systems  with  random  pinning,  the  motion  is  nearly  always 
plastic  at  driving  forces  just  above  the  critical  force;  that  is, 
the  vortices  move  at  different  velocities  in  different  parts  of 
the  system  [25].  Numerical  simulations  show  that  the  plastic 
motion  can  take  the  form  of  rivers  of  moving  vortices  sliding 
between  more  stationary  “river  banks”  [25]. 

4  Conclusion 

In  conclusion,  we  imaged  for  the  first  time  the 
phase  evolution  of  steady  state  vortex  motion  driven  by  an  ac 
transport  current.  A  quantitative  study  by  TRMO  imaging  and 
FEM  simulations  reveals  the  formation  and  coexistence,  dur¬ 
ing  the  ac  cycle,  of  a  quasi-static  glassy  state  in  the  sample 


interior  and  a  disordered  dynamic  state  near  the  edges  of  the 
sample. 
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